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ABSTRACT 

Star forming regions are expected to show linear proper motions due to the relative motion of the Sun with 
respect to the region. These proper motions appear superposed to the proper motions expected in features 
associated with mass ejection from the young stellar objects embedded in them. Therefore, it is necessary to 
have a good knowledge of the proper motions of the region as a whole in order to correctly interpret the motions 
associated with mass ejection. In this paper we present the first direct measurement of proper motions of the 
NGC 1333 star forming region. This region harbors one of the most studied Herbig-Haro systems, HH 7-11, 
whose exciting source remains unclear Using VLA A configuration data at 3.6 cm taken over 10 years, we 
have been able to measure the absolute proper motions of four thermal sources embedded in NGC 1333. From 
our results we have derived the mean proper motions of the NGC 1333 star forming region to be /i^ cos 6 = 9 
± 1 mas yr"' and fis = -10 ± 2 mas yr"' . In this paper, we also discuss the possible implications of our results 
in the identification of the outflow exciting sources. 

Subject headings: astrometry — ISM: individual (NGC 1333) — ISM: jets and outflows — radio continuum: 
ISM — stars: formation — stars: individual (SVS 13) 



1. INTRODUCTION 

The stars in nearby regions of star formation are expected 
to show linear proper motions in the order of milliarcseconds 
(mas) per year. These proper motions are the result of the 
relative motion between the Sun and the stars studied, and can 
be obtained, to first order, by comparing the observed position 
of the star with respect to the reference frame of the remote 
quasars at different epochs. 

In the last few years. Very Large Array (VLA) observations 
taken with time baselines of up to two decades have been 
used to measure these proper motions for regions like Tau- 
rus (Loinard et al. 2003), L1527 (Loinard et al. 2002), L1551 
(Rodriguez et al. 2003), Ophiuchus (Curiel et al. 2003), and 
Orion (Gomez et al. 2005). 

Even when the VLA observations lack sufficient angular 
resolution to detect more subtle motions (i.e. geometric par- 
allax), they can provide not only the average proper motion of 
the region but also reveal the presence of stars moving with 
peculiar velocities with respect to the region, such as those 
found in the Orion BN/KL region (Rodriguez et al. 2005, 
Gomez et al. 2005), that suggest a runaway nature for some of 
the sources. Much more accurate astrometry can be achieved 
using Very Long Baseline Interferometry (VLBI) techniques 
(e.g. Loinard etal. 2007), but these observations are restricted 
to very compact and bright non thermal radio stars and cannot 
be applied to thermal sources, where the emission is relatively 
extended and its brightness temperature is not expected to ex- 
ceed 10"^ K, far below the sensitivity of present day VLBI. 

The NGC 1333 star forming region, located at a distance of 
235 ± 18 pc (Cernis 1990; Hirota et al. 2008) in the Perseus 
complex, harbors the classical bright Herbig-Haro (HH) sys- 
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tem HH 7-11, first reported by Herbig (1974) and by Strom et 
al. (1974). The optically visible star SVS 13, discovered as 
a near-IR source (Strom et al. 1976), is roughly aligned with 
the chain of HH objects and was proposed as the powering 
source of this HH system. This association was questioned 
by Rodriguez et al. (1997), who discovered a cm radio source 
(VLA 3) located 6" to the SW of SVS 13, and ai-gued that 
this new object is a better candidate to drive the HH outflow. 
SVS 13 is also associated with cm emission (source VLA 4 
of Rodriguez et al. 1997) and mm emission (Looney et al. 
2000). Bachiller et al. (2000), through interferometric obser- 
vations, found that SVS 13 is associated with an "extremely" 
high velocity molecular outflow, although "standard" velocity 
gas was found in the vicinity of both SVS 13 and VLA 3. Sub- 
arcsecond VLA observations at 3.6 cm and 7 mm by Anglada 
et al. (2000, 2004) revealed that the radio source associated 
with SVS 13 is actually a close binary with two components, 
VLA 4A and 4B, separated by 0"3 ('-65 AU). Interestingly, a 
detailed analysis of the positions and spectral energy distribu- 
tion in the cm-mm range suggests that the two sources have 
very different properties: VLA 4B appears to be associated 
with the observed strong mm emission, probably arising from 
a circumstellar dust disk, while VLA 4A is the counterpart of 
the visible star SVS 13 with dust emission absent or much less 
significant in this component. 

Since high velocity water masers are indicators of outflow 
activity and can be observed with high angular resolution, Ro- 
driguez et al. (2002) analyzed archive VLA data of the water 
masers associated with SVS 13 in order to investigate which 
of the two stars of the binary system was the most likely can- 
didate to drive the outflow. It was found that the water masers 
appear segregated in two groups, according to their position 
and LSR velocity. A group with the LSR velocity similar to 
that of the ambient cloud is associated with VLA 4A, and a 
blueshifted velocity group is associated with VLA 4B. This 
result was interpreted as favoring VLA 4B as the outflow driv- 
ing source. However, the study of Rodriguez et al. (2002) 
was restricted to the line-of-sight component of the velocity. 
Recently, Hirota et al. (2008) presented VLBI astrometric 
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TABLE 1 
Observations Parameters 









Bootstrapped 












VLA 


Flux Density of 


Synthesized Beam'' 








Observation 


Project 


Phase Calibrator" 


HPBW 


RA. 


rms Noise'' 


Epoch 


Date 


Code 


(Jy) 


(arcsec) 


(deg) 


{fiJy beam"' ) 


1989.1 


89-Jan-14 


AR202 


1.343 ± 0.003 


0.29 X 0.24 


62 


28 


1994.3 


94-Apr-23 


AR277 


1.686 ±0.006 


0.28 X 0.20 


-68 


27 


1996.9 


96-Dec-22 


AR277 


1.305 ± 0.009 


0.35 X 0.30 


-88 


24 


1998.2 


98-Mar-27 


AA218 


1.476 ± 0.009 


0.26 X 0.24 


-5 


18 


1998.4 


98-May-26 


AA218 


1.610 ±0.010 


0.26 X 0.24 


-7 


18 


1999.5 


99-Jul-03 


AA239 


1.681 ± 0.008 


0.29 X 0.25 


39 


17 



" The phase calibrator used in all the epochs was 0333+321. From naturally weighted maps. 



TABLE 2 
Positions of the Sources" 



VLA 2 



VLA 3 



Epoch 



1989.1 
1994.3 
1996.9 
1998.2 
1998.4 
1999.5 



a(J2000) 
03 29 01.9535 ±0.0003 
01.9586 ±0.0004 

01.9605 ±0.0004 

01.9606 ±0.0002 
01.9604 ±0.0002 
01.9622 ±0.0002 



i(J2000) 



31 15 38.307 ± 0.005 
38.221 ± 0.004 
38.285 ± 0.008 
38.188 ± 0.005 
38.211 ±0.006 
38.183 ± 0.003 



a(J2000) 
03 29 03.369 ± 0.002 

03.375 ± 0.002 

03.376 ± 0.001 
03.376 ± 0.001 
03.373 ± 0.002 
03.373 ± 0.001 



6(nm)) 

31 16 01.89 ±0.02 
01.82 ±0.02 
01.75 ±0.02 
01.72 ±0.02 
01.78 ±0.02 
01.78 ±0.01 



VLA 4A VLA 4B 



Epoch 


Q(J2000) 


i(J2000) 


Q(J2000) 


i(J2000) 


1989.1 


03 29 03.730 ± 0.002 


31 16 04.02 ±0.02 


03 29 03.751 ±0.002 


31 16 04.08 ±0.02 


1994.3 


03.730 ± 0.002 


03.95 ± 0.02 


03.758 ± 0.004 


04.01 ± 0.02 


1996.9 


03.731 ±0.002 


03.94 ± 0.02 


03.757 ± 0.002 


04.01 ± 0.02 


1998.2 


. . .b 


. . .b 


. . .b 


. . .b 


1998.4 


03.735 ± 0.002 


04.00 ± 0.03 


03.760 ± 0.002 


03.99 ± 0.02 


1999.5 


03.734 ± 0.002 


03.96 ± 0.02 


03.757 ± 0.001 


03.91 ± 0.02 



" Positions derived from Gaussian ellipsoid fits in the naturally weighted maps. Units of right ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes, 
and arcseconds.'' The components of the binary system VLA 4 could not be resolved in this epoch. 



observations with VERA and derived absolute positions and 
proper motions of the order of 15-20 mas yr"' (correspond- 
ing to velocities of 16-22 km s"') of the water masers associ- 
ated with SVS 13. These observations confirm the positional 
association with VLA 4A of the masers with an LSR veloc- 
ity close to that of the ambient cloud. However, the proper 
motions (velocity and direction) relative to the sources could 
not be firmly established because of a poor knowledge of the 
absolute proper motions of the region. In fact, the optical 
data of Herbig & Jones (1983) suggest that the stars pro- 
jected upon heavy obscuration (presumably associated with 
the NGC 1333 cloud) have an average proper motion of ^10 
mas yr"' relative to the stars near the edge of the obscura- 
tion (presumably background stars). This suggests that the 
absolute proper motion of the region could have a significant 
contribution to the proper motions of the masers observed by 
Hirota et al. (2008). 

In this paper we present an absolute astrometry analysis and 
proper motion calculations over 10 years of four thermal ra- 
dio sources in the NGC 1333 region. There are no reports 
of absolute proper motions for embedded young stars in this 
region and our study provides the first determination of these 
motions. 



Data at 3.6 cm continuum were taken using the VLA of 
the National Radio Astronomy Observatory (NRAO)'* . The 
data consist of six epochs of observations, ranging from 1989 
January 14 to 1999 July 3, obtained using the A configu- 
ration of the VLA, which provides an angular resolution of 
^0"3. A requirement to obtain reliable astrometry is that the 
same phase calibrator has been used for most or preferably all 
observations. This requirement is accomplished in our case 
since all observations were made using 0333H-321 as phase 
calibrator. Data editing and calibration were carried out using 
the Astronomical Image Processing System (AIPS) package 
of NRAO, following the standard VLA procedures. The ob- 
servation dates, the bootstrapped flux densities of 0333H-321, 
as well as the parameters of the synthesized beams, and the 
rms noise of the naturally weighted maps are given in Table[T] 
Since the nominal position of the phase calibrator used in the 
first epoch was slightly different from that used in the other 
five epochs (the positions of the phase calibrators are period- 
ically refined at the VLA), in order to obtain more accurate 
astrometry, we corrected its position using the task CLCOR 
of AIPS before calibration of the first epoch. Positions of the 
four thermal sources (VLA 2, VLA 3, VLA 4A, and VLA 4B; 
see Rodriguez et al. 1997 and Anglada et al. 2000 for the 
nomenclature of the sources) were derived in each epoch by 
Gaussian ellipsoid fits in the naturally weighted maps. 



2. OBSERVATIONS 



'' The NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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Fig. 1 . — Naturally weighted VLA 3.6 cm continuum map of the thermal sources VLA 2, VLA 3, VLA 4A and VLA 4B at the 1998.4 epoch. Contours ai'e -3, 
3, 4, 5, 6, 8, 10, 12, 14, 16, 18, 20, 25 and 30 times the rms of the map, 18 fijy beam"'. The synthesized beam, shown in the right hand panel, is Of'26 X Of'24; 
P. A. =—7°. The three panels are plotted at the same scale. 



TABLE 3 
Proper Motions of the sources 





fict cos 6 






V 


PA. 


Source 


(mas yr"' ) 


(mas yr"' ) 


(mas yr"') 


(kms-') 


(deg) 


VLA 2 


10 ± 1 


-10 ±4 


14 ±3 


16 ±3 


140 ± 10 


VLA 3 


6±3 


-13 ±4 


14 ±4 


16 ±4 


160 ± 10 


VLA4A 


6±3 


-4±4 


7±3 


8±3 


120 ± 30 


VLA4B 


9±3 


-13 ±4 


16 ±4 


18±4 


150 ± 10 


Weighted Mean" 


9± 1 


-10 ± 2 


14 ±2 


16 ±2 


140 ± 10 





" The values of the total proper motion, velocity in the plane of the sky, and P.A. are derived from the mean values of /ia COS 5 and fis ■ 



3. RESULTS AND DISCUSSION 

In Figure[T]we show the 1998.4 maps of the sources VLA 2, 
VLA 3, VLA 4A, and VLA 4B. The measured absolute po- 
sitions of each source at each epoch are given in Table |2] In 
Figure|2]we present position vs time diagrams for each source, 
and the proper motions, fi^cosS and /i^, obtained by linear 
least squares fits of these diagrams for each source are given 
in Table[3] All the sources show proper motions roughly to the 
SE with P.A. in the range from --120° to --160° and veloci- 
ties in the range from 16 to 18 mas yr~', except for VLA 4A, 
whose velocity is about one half the velocity obtained for 
the other sources. This possible velocity difference of 8-10 



km s ' between VLA 4A and the ambient cloud may account 
for the yet unexplained large increase in the optical and near- 
IR brightness of VLA 4A that took place in 1988-1990 (Eis- 
loffel et al. 1991). The source may have simply displaced 
from a very obscured region to a position where it became 
more easily detectable. However, because of the signal-to- 
noise limitation of the data plus the errors produced by the in- 
trinsic differences in the morphologies, orientations and flux 
densities of the sources (see Fig. 1), we cannot conclude that 
the differences between the values obtained for each source 
indicate statistically significant differences in the motions of 
the objects. 
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Fig. 2. — Position vs time diagrams of the sources VLA 2, VLA 3, VLA 4A, and VLA 4B. The positions of each source are represented in the diagrams as 
right ascension and declination offsets (in arcsecs) relative to the mean of the positions given in Table|2] En'or bars are obtained from Gaussian ellipsoid fits in 
the naturally weighted maps. The actual en'ors could be greater because the moiphology of the sources is not exactly Gaussian (see text for details). The dashed 
line in each panel is a linear least squares fit to the data. The values of /la cos 5 and fig obtained from these fits are also labeled in each panel. 



As a representative value of the proper motion of the 
NGC 1333 region as a whole we adopt the weighted aver- 
age of the proper motion of the four sources, /i^ cos 5 = 9 ± I 
mas yr~' and = -10 ± 2 mas yr"'. Assuming a distance of 
235 pc for NGC 1333, these values imply a velocity of 16 ± 2 
km s~' along a P.A. of 140° ± 10°. This is the first measure- 
ment of the proper motion directly made on embedded objects 
that belong to the NGC 1333 star forming region. Therefore, 
we believe that our value is, up to date, the most representative 



of the proper motion of NGC 1333. 

The proper motion of the NGC 1333 region obtained from 
our data can be compared with the value expected due to the 
effect of the differential galactic rotation and the motion of 
the Sun with respect to the LSR. Assuming the solar mo- 
tion with respect to the LSR based on the Hipparcos satellite, 
(Uo, Vo, Wo)=(10.00, 5.25, 7.17) km s"' (Dehnen & Binney 
1998), the model of Brand & BHtz (1993) for the galactic ro- 
tation, and a distance of 235 pc for NGC 1333, we would ex- 
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pect proper motions of /i^ cos 5 = 2.9 mas yr~' and fis = -7.5 
mas yr~', that imply a velocity of ^9 km s"' along a P. A. of 
~160°. This velocity is ~7 km s~' smaller than the value ob- 
tained from our data. Since the difference between the proper 
motion predicted by the model and our measurements is of 
the order of the one-dimensional rms velocity dispersion ob- 
served between molecular clouds (7.8 km s~'; Stark & Brand 
1989), we attribute this difference to the peculiar motion of 
the NGC 1333 cloud. We have looked for optical stars as- 
sociated with NGC 1333 with proper motions measured by 
Hipparcos. The only nearby star we find with these character- 
istics is BD + 30°549, a B8:p type star that illuminates part of 
the reflection nebula NGC 1333 (e.g. Rodriguez et al. 1990). 
For this star, Hipparcos reports proper motions of /Zq, cos 5 = 
6.31 ±2.35 mas yr~^ and ns = -8.87±1.77 mas yr~^ (Perryman 
et al. 1997). These values agree with the weighted average 
obtained by us at the ±l-a level. 

Our results confirm that the source VLA 3 is a member of 
the NGC 1333 molecular cloud since this source shows the 
same proper motion that the other members of the cloud. Al- 
though this conclusion alone does not resolve the problem of 
the identification of the powering source of the HH 7-11 sys- 
tem, our results discard the possibility that VLA 3 is a fore- 
ground or background source not related to the NGC 1333 
region. Therefore, the proximity of VLA 3 to the base of the 
HH 7-11 system is physically real, and not only a projection 
effect. 

The determination of the proper motions of embedded 
young stars in NGC 1333 is useful to reanalyze the proper 
motions of water masers reported by Hirota et al. (2008). 
These authors report VLBI observations of two maser fea- 
tures associated with VLA 4A (SVS 13). They found absolute 
proper motions of the order of 15-20 mas yr~' in a direction 
(P.A.~114°-133°) roughly ahgned with the jets and outflows 
in the region. However, as noted by Hirota et al. (2008), these 
proper motions contain a significant contribution due to the 
relative motion of the NGC 1333 region with respect to the 
Sun. In order to correct for the contribution due to the rela- 
tive motion of the NGC 1333 region with respect to the Sun, 
Hirota et al. (2008) subttacted the theoretical proper motion 
of the region due to the Solar motion relative to the LSR (ne- 
glecting the contribution of Galactic rotation, that they show 
is small). This correction results in proper motions of the or- 
der of 7-14 mas yr~', and a significant change in the direction. 
The proper motions corrected in this way are not parallel to 
the HH 7-11 jet and CO outflow. However, as we have dis- 
cussed above, there is still a significant contribution due to 
the peculiar proper motion of the NGC 1333 region that must 
be taken into account for a proper interpretation of the results. 

Our results show that the mean proper motion of the 
NGC 1333 region (/i„cos5 = 9±1 mas yr~^ and /x^ = 
— 10±2 mas yr ) determined from the four embedded thermal 
sources is comparable to the absolute proper motions found 
by Hirota et al. (2008) for the two water maser features as- 
sociated with VLA 4A (/Zc cos 6 = 17.9±0.9 mas yr"' and 
= -7.9±1.4 mas yr~' for feature 1, and /ic«cos5 = 10.6±1.7 
mas yr~' and = -10.0±2.1 mas yr~' for feature 2). There- 
fore, this correction appears to be very important to obtain the 
true proper motions relative to the region. The close coinci- 
dence for feature 2 suggests that this maser is almost station- 
ary with respect to the cloud and that is not ttacing sigiuficant 
peculiar proper motions. This conclusion is supported by the 
fact that the LSR velocity of the masers (7-8 km s~') nearly 
coincides with the LSR velocity of the molecular cloud, as 



already noted by Rodriguez et al. (2002). After correction 
for the proper motion of the NGC 1333 region determined 
by us, feature 1 has residual proper motions of ^aCos5 = 
9±2 mas yr~' and fis = 2±3 mas yr"^ This indicate motions 
to the east (P.A.~90°) wifli a velocity of ~10 km s~\ that 
align poorly with the chain of bright HH 7-11 objects and the 
molecular outflow. In any case, these small residual proper 
motions align much better with the faint chain of HH objects 
(a, b, c and d) found by Davis et al. (1995), that aligns in the 
west-east direction (see Fig. 2 of Rodriguez et al. 1997 where 
the positions of the sources discussed in this paper and the HH 
objects are also shown). 

We note, from the analysis discussed above, that we ob- 
tained proper motions of the maser features relative to the 
NGC1333 cloud. An unambiguous measurement of the 
proper motions relative to VLA 4A would require to subtract 
the individual proper motion of this source instead of the av- 
erage motion of the NGC 1333 region. The nominal proper 
motion velocity we obtained for this source is about one half 
the average value of the region. Taken at face value, it could 
suggest that the VLA 4A source has a significant velocity rel- 
ative to the ambient cloud. However, the weakness of the 
source does not allow a reliable determination of its individual 
proper motion with the currently available data. Additional, 
high sensitivity data would be required to investigate possi- 
ble motions of the VLA 4A source relative to the NGC 1333 
cloud in order to accurately determine the motions of the wa- 
ter maser features relative to VLA 4A. 

4. CONCLUSIONS 

We have presented the results of multi-epoch high angu- 
lar resolution VLA observations at 3.6 cm of the NGC 1333 
star forming region. From these observations we have mea- 
sured the absolute positions at different epochs of four thermal 
sources (VLA 2, VLA 3, VLA 4A and VLA 4B) tracing YSOs 
embedded in this region. All these sources show an absolute 
proper motion to the SB, which we interpret as a proper mo- 
tion of the NGC 1333 region as a whole. From our data we 
obtained the average values fiaCosS = 9 ± 1 mas yr"' and 
Us = -10 ± 2 mas yr"\ that we interpret as representative of 
the proper motion of the molecular cloud associated with the 
NGC 1333 region as a whole. 

Our results allow us to correct the proper motion measure- 
ments of water maser spots associated with VLA 4A (SVS 13) 
obtained by Hirota et al. (2008) in order to obtain the proper 
motions relative to the NGC 1333 cloud. We conclude that 
feature 2 of Hirota et al. (2008) is practically stationary with 
respect to the cloud, and that the small residual proper mo- 
tions of feature 1 indicate motions to the east, suggesting that 
none of these features is associated with the HH 7-11 system 
or the molecular outflow. 

Our data also suggest that VLA 3 is a member of the 
NGC 1333 molecular cloud, since this source shows a proper 
motion similar to that of the others sources embedded in the 
region. Since the position of VLA 3 appears close to the base 
of the HH 7-11 system, we conclude that this proximity is 
physically real, and not a projection effect. 
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